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Abstract 



H 

^' 

jrt ■ Strong mixing property holds for a broad class of linear and nonlinear time series 

models such as ARM A and GARCH models. In this article we study correlation 

structure of strong mixing sequences, and some asymptotic properties are presented. 

^ . We also present a new method for detecting change point in correlation structure 

Q^ ■ of strong mixing sequences, and present a nonparametric CUSUM test statistic 

f— ^ I for this. Asymptotic consistency of this test statistics is shown. This method is 

cn ! applied to simulated data of some linear and nonlinear models and power of the 

O ■ 

p^ ■ test is evaluated. For linear models, it is shown that this method have a better 

performance in compare to Berkes et al.(2009). 
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1 Introduction 

Change point detection in a sequence of random variables was first proposed by 
Page(1954). This study started by detecting changes in the mean of a sequence of in- 
dependent random variables and then extended to dependent sequences. Change point 
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detection is widely used in various fields such as quality control, economics, finance and 
medicine. Review of earlier works can be found in Csorgo and Horvath(1988), Brodsky 
and Darkhovsky (1993) and Csorgo and Horvath(1997). 

Among different methods for change point detection, the CUSUM test proposed by 
Page(1954), for mean change detection, is widely used for its simplicity. Inclan and 
Tiao(1994) proposed a CUSUM of squares test for testing a variance change in i.i.d. 
normal random variables. Lee and Park(2001) extended the CUSUM test of squares 
test of Inclan and Tiao(1994) for linear processes. Lee et al.(2003) studied change of 
parameters in a random coefficient AR(1) model, thus detecting changes in the auto- 
covariances of a linear process Galeano and Pena(2007) studied changes in variance and 
correlation structure of the multivariate time series. Zhou and Liu(2009) used a weighted 
CUSUM statistic for mean change detection in infinite variance AR(p) process. Berkes et 
al.(2009) considered a CUSUM test to detect changes in the mean and in the covariance 
structure of a linear process. Recently Qin et al.(2010) studied mean change detection in 
a-mixing processes. 

In this article we study change in the correlation structure of strong mixing sequences. 
Let {Xt : t > 1} be a stationary process. As a measure of dependence we use Rosenblatt's 
a-mixing coefficient as 

ax{n,j) = sup \PiA n B) - P(A)P(5)|, (1.1) 

A,B 

where A and B are in the a-fields M^^(X) = (T{Xt;t < n} and M^_^^(X) = cr{Xt;t > 
n + j} respectively. The sequence {Xt} is said to be a-mixing or strong mixing (SM) if 

axij) = supaxin,j) ^0 as j ^ oo. (1.2) 

n 

Strong mixing processes are asymptotically independent. Strong mixing property holds 
for a large class of linear and nonlinear stationary time series such as ARMA and GARCH 
models, m-dependent processes, broad class of Gaussian processes and ergodic Markov 
processes( Bosq 1996, and Bradley 2005). 

Ibragimov(1962) showed some results for stationary strong mixing sequences and 
proved central limit theorem for strict stationary SM processes. Davydov(1968) ob- 
tained some moment inequalities and Rio(1993) presented some covariance inequalities 
and bounds on the variance of partial sums of SM processes. IIerrndorf(1985), Doukhan 
et al.(1994), and Merlevede and Peligard(2000) studied functional central limit theorem 



on SM processes. Romano and Thonibs(1996) used central limit theorem, established 
by Ibragimov(1962), to show that sample auto-covariances of strictly stationary SM se- 
quences converge in distribution to normal distribution. 

By using functional central limit theorem for SM sequences, we propose a new test 
statistic for detecting changes in correlation structure of stationary strong mixing pro- 
cesses. This test is a nonparametric one and does not depend on any assumptions about 
the underlying distribution or model. The rest of this paper is organized as follows: In 
section 2 a nonparametric test statistic for detection of change points, in broad class 
of linear and non linear process, is constructed and its asymptotic properties, under no 
change null hypothesis, are studied. Also the consistency of this test statistic is shown in 
section 2. Section 3 is devoted to the simulation results on different linear and nonlinear 
models. In this section the method of Berkes et al.(2009) has been compared with the 
method of this paper for some some linear models. By simulation we show that this test 
statistic have a better performance and is more powerful in many cases. 

2 Main results 

In this section we present some preliminary results which will be used later in this 
paper. We present functional central limit theorem for sample auto-covariances of SM 
processes. We also introduce a new test statistic for detecting changes in the correla- 
tion structure of stationary SM processes, which we call CUSUM strong mixing(CSSM). 
Finally we show consistency and asymptotic convergence of this empirical CSSM test 
statistic. 

Let {Xn} be a sequence of random variables on some probability space {Q,J^,V), 
satisfying 

E(X„) = 0, E{Xl) < oo for n > 0. (2.1) 

Let Sn = Xi + ... + Xn for n > 0. Consider the Skorokhod space D = D[0, 1] of all 

functions on [0, 1] which are right continuous with left limit. Let Wn{t) : Q ^ D to he a. 

random function as 

WJt) = -^ for te [0,1], n>0. 
a^/n 

If Wn{t) is weakly convergent to a standard Brownian motion W{t), then X„ is said 

to satisfy the functional central limit theorem or strong invariance principle (Billings- 



ley, 1999). Herrndorf(1985) proved functional central limit theorem for strong mixing 
sequences without stationarity assumption but assumed convergence of the variance of 
partial sums. 

For Xi, X2, ... as a sequence of zero-mean stationary process the sample auto-covariances 
7„(/i), h = 0, ...,n, are defined as: 

1 n—h 

ln{h) = 7„(-/i) = - E X^X,+h. (2.2) 

Asymptotic covariance of sample autocovariances is known as Bartlett's estimator and is 
defined as: 

Chk = lim^nCov{%{h),jn{k)). (2.3) 

Let 

g:{t) := \^C-'/\^^^,^{0) - 7n(0), 7h(1) - %{!), -, llnt]{L) - 7„(L))^, (2.4) 



where C = [chk]hX=i is the covariance matrix whose entries are defined by (2.3). Now we 

have the following result. 

Theorem 1: Let {X„}5^^ he a stationary strong mixing process which satisfies: 

• z;supiE|Xi|4+2'5<oo 

• n) ET=i Cix{ky'^'^+^'^ < 00 , for some 6 E (0, 00) 

in which ax{-) is the mixing coefficient, defined by (1.2). Then 



L 

(\A 



9*n{tf -at) ^ E(W^°(^))' asn^^, (2.5) 

i=o 



where =^ denotes convergence in distribution, Wj'(.) are independent Brownian bridge, 
for 0<j<L. 

Before proceeding to the proof of this theorem we present some lemmas, which are nec- 
essary for our proof. 

Lemma 1: (Davydov, 1968) Let the process {Xt} be strong mixing, and random vari- 
ables C, and rj be measurable with respect to M"q^(X) and M^(X), introduced by (1-1), 



respectively. Moreover if the moments E\C^\^ and E\ri\'^ exist for p,q > where - + - < 1, 
then 

\ECv - ECEr]\ < C[E\^\PY/P[E\r]\^Y/''[a{n)Y-^/P-^/''. 

Lemma 2: Let {Xn}'^=i be a zero-mean strong mixing process, where 

oo 

^'5/(2+'5) ^ ^ „^,j ^„^ f^\ ^.|2+5 = M <oo, for some 5 E (0, oo). 



J2c^x{ky'^'^^'<oo, and snp E\Xi 
fc=i « 



Then 



EjSr,? 



is convergent. 



Proof of Lemma 2: By lemma 1, 

t=l s=l t=l s=l 

Let p = q = 2 + 5, then 

-| 1 n CO 

-E{S^f < -M^ E n[axit)f'^' < M E [^xit)f'''' < oo. 

i=—n i=—oo 

An alternative proof for lemma 2 can be found in Rio (1993). 

Proof of Theorem 1: Let 

grXt) := ^C-^/2(7h(0) -7(0),7m(1) -7(1),-,7h(1^) -7(i^))^, (2.6) 

where 7(.) is auto-covariance function of {Xf}. By (2.4) and (2.6), it is immediate that 

g:{t)=gn{t)-^^gn{l). 
n 

Any Brownian bridge W^{t) has the same distribution as W{t) — tl^(l), where W{t) is 

a standard Brownian motion (Billingsley, 1999). 

Therefore ii gnit) =^ (Wo{t),Wi(t), ...,WL(t))^, where Wi(t) are independent Brownian 

motions for < i < L, then g^t) =^ (W^it), W^{t), ..., Wl{t)f. So it is enough to show 

that gn{t) converges to a vector of Brownian motions. 

As 



/ 



gn{t) 



M^-1 



^-1/2 



n 



7[nt](0)-7(0) 



\ 



C-V2 



\ 7[„t](L)-7(L) J 



n 



( eCU.^. - K]7(o) ^ 



VEtt^x,x,+^-K]7(L); 



so by assuming 



we have that 



Ym,t = XtXt+rn - lim), for <m < L 



(2.7) 



9nit) 



^-1/2 



n 



( Y^H Y \ 






^-1/2 



n 



i \ 

7(1) 
27(2) 

V ^7(^) ] 



{2.1 



By (2.7), 



and 



^-ooiY) = ^{Ym,u t<n} = a{X,X,+„; t < n} C MI^^Tl^) 



M~ ^.(F) = a{Yrn,u t>n + j} = a{X,X,+„; t > n + j} C M^ ^.(X). 

So relation (1.1) imphes that, ay^ij) < ax(j — "o^)- Hence {Vm,t} form a zero mean, 
strong mixing process, where by assumption (ii). 



E«yJA:)V(2+^)<oo. 



k=l 



Also by (2.7) and assumption (i) 



2+5 



sup E\Ym,i\ < OO 



(2.9) 



(2.10) 



for some 6 G (0, oo). Let Sm,n = ^m,i + ••• + ym,n- Using (2.9) and (2.10), lemma 2 implies 
that ^°''^f"-"^ ^ o-^, for some o"^ < oo. 

If (Tm > 0, then (2.9), (2.10), and functional central limit theorem, introduced by Herrn- 
dorf(1985), assert that 

'-'m,[nt] 



crm.\/n 



FoTO<h,k< L, 



[nt]-h 



[nt]-k 



W{t), for 0<m<L. 



[nt]-h 



(2.11) 



[nt]-k 



cov 



— J2 Yh,^,^ E >^m) = cov(— J2 Yh,+^h^{h),^ J2 Yk,,+^k^{k)) 



1=1 



^ 1=1 



^ 1=1 



n 



n 



i=l 



n 



nt]'^ \nt] 

cov(7[„i](/i) -7(/i),7[„i](A;) - 7(A;)) = [nt]coY{%t]{h),^[nt]{k)). 



n ■ ■ ■ n 

By Bartlett's formula (2.3), hm„_^oo[nt]cov(7[„(](/i),7[„(](A;)) = Chk- So 



1 



[nt]-h 



1 



[nt]-k 



lim cov(^ V Yh,i,—;= V Yk,i)=tChk- 



i=l 



^ ;^i 



(2.12) 



Hence by (2.11) and (2.12), 

^-1/2 



^ EHV™. \ ( Wo{t) 



n 



1 Y^["t]--^v 

y 2^1=1 ^m,i ) 



\ W^L(t) ) 



(2.13) 



where Wj{t\2 = 0,...,L, are independent Brownian motions. 

The second part on right hand of (2.8) tends to zero as n — )■ oo, so (2.13) imphes that 



/ 



9n{t) 



Woit) 



\ 



V WLit) ) 



where Wj^j = 0, ...,L, are independent Brownian motions. □ 

2.1 CUSUM test statistic 

Using theorem 1 and (2.4), a CUSUM test statistic is constructed as: 



T„ := max q 

L<k<n '--n 



:{-Y-9V 



n 



(2.14) 



By continuous mapping theorem 



Q(+\\'i 



T„ ^ sup E(W^°W) 



0<t<l 



(2.15) 



i=o 



For detecting changes in time series {Xt}, under the assumptions of theorem 1, the 
following test is proposed for testing hypothesis 

Hq. no change occur in the auto-covariance function of Xi, ...,X„ 

Hi. there is a 1 < A; < n such that auto-covariance function of Xi, ..., Xk is different from 
auto-covariance function of X^+i, ..., X„. 

The strategy of this test is to reject Hq when T„ is large. 

By (2.14) and (2.15), the critical region of the test at significant level a is {T„ > Cq}, 
where c^ is the (1 — a)-quantile point of the distribution of suPq<j<]^ Z]j=o(^j°(^))^- The 
critical values can be found in Kiefer(1959) and Lee et al.(2003). 

Example: Let {Xt} be an MA(1) process as: 



Xt = Zt + 9Zt-i, 



r igUre li behavior of T„ for different values of parameter 9 in MA(1) model 




where {Zt} is an iid normal sequence with mean zero and variance a^. In hnear processes 
where E[Zf) = rja'^, Bartlett's formula has explicit form as: 

where 7(/) is the corresponding autocovariance function at lag / of {^t}, see Brockwell 

and Davis(1991). 

If the noise is Gaussian, 77 = 3, so 

Ci,i = 7(l)7(j - i - 1) + 7(1 + i)l{3 - 1) + 7(0)7(j - i) + li^hU) + li^hiJ - « + 1) 

+7(l+j)7(l-^)- 
Let L = 1, therefore corresponding covariance matrix C = [cij]f j^i can be written as 

272(0) +472(1) 47(0)7(1) \ / 2(1 + 4^2^^4) 4^(1 + ^2) 



C 

47(0)7(1) 7^(0) + 37^(1) 

By (2.4), for t = ^ we have that 

k 



40(1 + 02 



:i + 502 + e 



2 I a4\ 



a'. 



gl{k/n) := —C-'I\%{Q) - %{0), %{1) - %{l)y 



n 



and by (2.15) 



e 



T„ = max {-(7fe(0) - 7„(0), 7^(1) - %{l))C-\%{Q) - 7„(0), 7^(1) - 7n(l))^}. 

L<k<n n 

Figure 1 shows the behavior of test statistic T„ for different values of parameter 6 in 
an MA(1) process without change point. For L = 1, at significant level a = 5%, the 
critical value is c^ = 2.408. 
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2.2 Consistent estimation of covariance matrix 

In this section we evaluate asymptotic behavior of the covariance function of estimators 
7„(i), defined by (2.3). Bartlett(1946) derived an exphcit formula for the asymptotic 
behavior of covariance function of sample autocovariances, when there exist a linear 
model for the data(Priestley, 1981). We present a consistent estimator C for the case 
that there is no model for data, or we have nonlinear process. Using (2.3), Bartlett's 
estimator can be written as: 

Chk = lim 9n(h,k) = lim ncov(%(h),%(k)), < h < k < n. 

n-^■oo ^ ' n-^oo ^ ^ ' ^ '' 

For stationary process {X^}, let 

1 " 
7n(/i) = -E^*^t+^. 0</i<n, (2.16) 

and 6n{h,k) = 'ricov(7„(/i),7„(A;)). By theorem 3 it is shown that 6n{h,k) has the same 
asymptotic behavior as 6n{h, k). For the evaluation of dn{h^ k), one can easily verify that 

9n{h,k) = ncov(7„(/i),7„(A;)) = n{E{%{h).%{k)) -j{h).-f{k)} 

-inn 
■inn 1 n— 1 n—l 

= -T.T. E{XtXt+HXsX,+k - i{hh{k)} = - E T^i^i^tXt+hXt+iXt+i+k) 

n — — 11 — — 

"' t=l s=l "' 1=0 t=l 

-I —1 n 

-7(/^)7(A;)} + - E E {EiXtXt+,XmXt+i+k) -lihhik)}. 

^ l=~n+lt=-l+l 

By replacing / with — / and t — / with t in last summation we have 

^~„(/i,A;) = -E^m(/), (2.17) 

^1=0 

where 

n 

^h,m = T.{E{XtXt+hXtXt+k) - lih)^ik)}, (2.18) 

t=i 

and for 1 < / < n — 1 

a,Al) = E{^(n*) + EiY^t) - ^hhik)}, (2.19) 

in which F/^ = XtXt+hXt+iXt+i+k and Y^ = Xt+iXt+i+hXtXt+u- 



Let 

Uh,k) = -Y,ai,4l), (2.20) 



^1=0 



where 

n 1 n 



^k,m = E{- E{n° - %ih)%ik)}, (2.21) 



t=r^»=i 



n— i 1 n 



^/.,;i(0 = E{- Ein' + >^2' } - 27„(/i)7n(A;)}, l<l<n (2.22) 

and {hn} is a sequence of positive integers that 

/i„ = 0{n^) for some /3 G (0, 1/2). (2.23) 

Now we have the following result. 

Theorem 2: Under the assumptions of theorem 1, if sup^-E|X(|^+'^ < oo for some 5 > 

then 

II Onih^ k) — 9n{h, A;) II2 — )> as n —)■ 00. 

Proof : The proof is organized in three steps: 

Step 1: lim„_^oo^Ei<«<n|o'/i,fc(OI < 00. 
By lemma 1 and (2.19), 

n—l n—l 

t=i t=i 

n—l 



n—l 



Y^[E\Xt^iXt+iU''Y^^[E\XtXt^k\'Y^'^[ail + h)]'-'/^-'/'^. 
t=\ 

So for p = g = 2 + 5, by assumption (ii) of theorem 1, 

\^Kkil)\ < M{n - l){ax{l - h)^ + axil + h)^}, 

and 

lim - E WhAl)\< 00. (2.24) 

l<l<n 

Step 2: ||a(/) -a(/)||2 is of order 0{n^/^). 
For 1 < / < n, by (2.19) and (2.22) 

fi — / fi — I '^ Ti "1 71 

lk.,.(/)-a.,.(/)|h = ||E{^n*+^>^2t-27W7(A:)}-E{-E^i+-E^/-27nW7n(^^ 

10 



n—l in 1 



''■ 1=1 

n—l 

< Ei^i + 02 + 203}, 



t=i ^ iTi n .^1 

n— i 



t=l 



where a^ = ||^ ELiI^' - ^^'tllh, a, = ||i Er=i{>^2' - ^^^itlb and a, = \\%{h)%{k) • 
'j{h)'j{k)\\2. As by the assumption of the theorem sup^ ElYt]"^^^^^ < 00, so by lemma 2 



1 " 



^ i=i 



and ai = 0(n^^/^). Similarly 02 = 0(n^^/^). Also 

03 = \\%{h)%{k) -7(^)7(^)112 = \\%{h)%{k) -%{h)'j{k) + %{h)'y{k) -7(^)7(^)112 

< ||7n(/i)||2||7n(A;) - 7(A;)||2 + \\ln{h) - -f{h)h\i{k)\. (2.25) 

Also by lemma 2, 

\\%ml = l|-EX,X,+,l|^ = ii?(EX,X,+,)2 < oo, (2.26) 



^ "i ^ t=i 



and 



in in 

l7n(A;)-7(A:)||^ = ||-E^*^*+^-7(A:)||^ < -^E{J2X,X,^,-^{k)f = 0{n-'). (2.27) 



^ t=i ^ t=i 



Thus by relations (2.25), (2.26) and (2.27) we have that 03 = 0{n-^/^). 

By similar method, one can easily verify that ||a'/i^fc(0) — o'/jfc(0)||2 = 0(n^/^). Therefore, 

\\crkAl)-(rh,km2<0{n'^^). (2.28) 

Step 3: Steps 1 and 2 are applied to prove the main result. 
By Minkowski inequality we have 

in -I hn 

WUKk) - e4h,k)h = \\-Y.^hAO - -T.^h,km2 < 



-I h„ 1 " 

E W^hAO - ^hAm2 + - E I^m(OI- (2.29) 



(=1 l=h 



So by (2.23) and (2.28) 



1 hn 1 



-T.\WhA^)-^hAOh <-hnO{n'^^)^0, as n^ 00. (2.30) 



^ITi n 
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As hn — ;■ c>o, (2.24) implies that 

J2 -|^/i,fc(OI ^0, asn^oo. (2.31) 

h„<l<n 

Finally by (2.30) and (2.31), we arrive at the assertion of the theorem. □ 

Theorem 3: Under the assumptions of theorem 2, we have that \\On{h, k) —6n{h, k)\\2 — )> 

0, as n — 7> oo. 
Proof : As 

\\en{h,k)-en{h,k)h< \\en{h,k)-en{h,k)h + \\Uh,k)-en{h,k)h, 

so by theorem 2 the second part on the right tends to zero, and for the first part, by (2.2) 
and (2.16) 

\\9n{h,k) -9n{h,k)\\2 = n||cov(7„(/i),7„(A;)) - cov(7„(/i),7„(A;))||2 

-1 n—hn—k n n 

= -II E E E{XtXt+,X,X,+k - lih^k)} - E E E{XtXt+f,X,X,^k - 7(/i)7(A;)}||2 

^ t=l s=l t=l s=l 

in n 

= - E E \\{XtXt^hXsX,^k - lihhik)}h ^ 0, as n ^ oo.D 

t=n—h+l s=n—k-\-l 

Corollary 1: Under the assumptions of theorem 2, by choosing C = [chk]L+ixL+i, where 
Chk = 9n{h,k) defined by (2.20), C* is a consistent estimator of covariance matrix C in 
relation (2.4). 

Corollary 2: Under the assumptions of theorem 2, by (2.4), (2.14), (2.15), and corollary 

1, we have 

f„ = max g:ik/nf.g:ik/n) => sup E«(^))'' (2-32) 

where 






g:^Xs) ■■= ^C'-'/'(7[h(0) - 7n(0),7M(l) - 7n(l), ...,7M(i^) - ln{L)f , (2.33) 



n 



in which C = [On{h, k)], and 9n{h, k) is defined by (2.20). 

3 Simulation results 

In this section, we investigate the performance of the proposed test statistic T„, by a 
simulation study. As this test statistic is to be applied for linear and nonlinear models, 
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Table 1: Empirical power of T„ for ARMA{1,1), where the initial parameter (9o,4>o) = (0.1, 0.2). 
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0.295 


0.874 
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0.989 
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0.734 


0.977 


0.994 


0.998 


0.7 




0.935 


0.995 


0.999 


1.000 



we consider simulations of such classes of time series. Test statistics are evaluated by 
using relations (2.32) and (2.33) with L = 1 and relation (2.20) with /i„ = n^'^. 

For creating change in the covariance structure of time series, the parameters are 
changed at the midpoint of the series. Empirical powers are evaluated, and for the case 
that there is no change in data, probability of type I error is reported. 

The critical value of the test statistic at level a = 0.05 is c^ = 2.408(Lee at el. 2003). 
Simulations are repeated 1000 times, for the following linear and nonlinear models, to 
evaluated empirical powers. 
Linear models: 

• Model 1: 

ARMA(1,1): X,-0X,_i = Z, + ^Z,_i, 

• Model 2: 

MA(2): Xt = Zt + diZt-i + ^2^t-2, 

where {Zt} is a sequence of iid normal random variables with mean zero and variance 1. 
Table 1 reports results of the test for simulated data from ARMA(1,1), where 250 
samples are generated with (^Oj^o) = (0.1,0.2), and then 250 more samples for different 
values of (^i,0i) as reported in the table. In table 1, empirical powers are evaluated for 
cases where one or both parameters have changed. In all cases high empirical powers 
shows ability of this test statistic. The value pointed by * is type / error, empirical level, 
which is slightly below 0.05. 

As a comparison of CSSM with the Berkes method for linear models, mode 2 with 
^ = ^1 = ^2 is considered and 250 samples are generated with ^ = at first stage and 
then 250 more samples for some alternative 9, and empirical powers are plotted in figure 
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Figure 2: Empirical power of CSSM in comparison with Berkes et al.(2009) in MA(2) model, Xt = Zt + BZt-i+BZt- 
parameter changes from to alternative 9. 
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r igUre O'. Behavior of CSSM and Berkes statistic in simulated samples from MA(2) model, without change point. 
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r igUre 4; simulated samples from 2-depcndent model with a change point in fc* = 500, variance of Zt change from 1 
to 1.26(left), Corresponding test statistics(right). 
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2. As empirical powers show, CSSM test has a better performance, with smaller type I 
error. 

To visualize this, we follow a simulation of model 2, where there is no change in time 
series data. So we generate 500 samples from MA(2) with (^1,^2) = (0.3,0.3) once and 
again with (^1,^2) = (0.2,0.7). Then we evaluate CSSM and Berkes test statistics, and 
plot them with corresponding critical values, 2.408 and 1.36 in figure 3. As figure 3 shows 
Berkes method indicate a change point by mistake where CSSM statistic is far beyond 
such miss detection. 
Nonlinear models: 

As nonlinear models, we consider foUowings: 

• Model 3: 

2-dependent: Xt = ZtZt-iZt-2 

• Model 4: 

GARCH(1,1): Xt = ht.Zt where h^t = 00 + aX^^ + /3h^t_^ 

in which {Zf} is a sequence of iid normal random variables with mean zero and variance 
one, 0"^ = 1. 

Figure 4(left) shows 1000 generated samples of a 2-dependent process, model 3, with a 
change point at k* = 500, where variance of {Zf} changes from 1 to 1.26. Figure 4(right) 
shows the behavior of CSSM and Berkes test statistics. Corresponding critical values are 
2.408 and 1.36 respectively which are presented by horizontal lines. Figure 4(right) shows 
that the supremums of both statistics exceed corresponding critical values and a change 
in process is detected, but CSSM statistics is more precise, as CSSM detects the change 
at t=512 and Berkes statistic detects it at t=750. 

Table 2 reports the empirical power of the CSSM, for 2-dependent model, model 3. In 
table 2(a) change of the variance of {Zf}, from o"^ = 1 to alternative values is proposed. 
In table 2(b) change of the mean of {Zf}, from /ig = to alternative values, has been 
considered. 

Table 3 shows the empirical power of the CSSM statistics in a GARCH(1,1) model. 
These simulations are done for different values of n. Parameters initial values are {u, a, /3) = 
(0.5,0.1,0.2), and empirical powers for different alternatives are reported. Simulations 
show that the powers has significant increase with sample size. 
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iabl6 L'. Empirical power in 2-dependent model. 
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4 Conclusion 

In this article a nonparametric CUSUM test statistic is proposed for detecting struc- 
tural changes in strong mixing time series. Under a sufficient condition this test statistic 
converges in distribution to the supremum of the sum of independent standard Brownian 
bridges. This method covers a broad class of linear and nonlinear time series such as 
ARMA and GARCH models, m-dependent models and many others. Beside the wide 
applications, simulation results shows that our test statistic in comparison with Berkes 
et al.(2009) has a better performance and is more powerful. 
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